Density of states of continuous and discrete spin models: a case study 
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v ■ A relation between 0(n) lattice spin models and Ising models denned on the same lattice was 

■ recently put forward [L. Casetti, C. Nardini, and R. Nerattini, Phys. Rev. Lett. 106, 057208 (2011)]. 

Such a relation, inspired by an energy landscape analysis, implies that the density of states of an 
, 0(n) spin model on a lattice can be effectively approximated, at least close to the phase transition, 

in terms of the density of states of an Ising model defined on the same lattice and with the same 
interactions. In the present paper we show that such a relation exactly holds, albeit in a slightly 
^f*j , modified form, in the special cases of the mean-field XY model and of the one-dimensional XY 

model. We also discuss the possible consequences of this result for the general case. 
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03 ■ I. INTRODUCTION 

Simple models play a prominent role in theoretical physics and especially in statistical mechanics. Moreover, 
finding links, relations or mappings — either exact or approximate — between different models often allows a deeper 
Q ■ understanding of the models themselves and of the physics they describe. 

In a recent paper [l[ a relation between the microcanonical densities of states of continuous and discrete spin models 
was conjectured. More precisely, it was suggested that the density of states of an 0(n) classical spin model on a lattice 
can be approximated in terms of the density of states of the corresponding Ising model, i.e., an Ising model denned 
on the same lattice and with the same interactions. Such a relation is suggested by an "energy landscape" approach 
[H to the microcanonical thermodynamics of these models, the key observation being that all the configurations of 
, an Ising model on a lattice are stationary points of an 0(n) model Hamiltonian defined on the same lattice with the 
■ same interactions, for any n. The relation between the densities of states can be written as 

W W(e)«o;W(e)flW( e ), (1) 

CM' 

7— I ' where e is the energy density, i.e., e = E/N where E is the total energy and N the number of spins, ui (n > is the 
(<— v density of states of the 0(n) model, oj^ 1 ' the density of states of the corresponding Ising model and gW is a function 
' representing the volume of a neighborhood of the Ising configuration in the phase space of the O(n) model. The 
t-H ■ function is unknown in general, but since it comes from local integrals over a neighborhood of the phase space 
one expects it is regular. The relation (fTJ) will be more precisely discussed in Sec. HTJ it is an approximate one and 
the approximations involved are not easily controlled in general. As discussed in PJ, were this relation exact there 
would be a very interesting consequence: the critical energy densities of the phase transitions of all the 0(n) models 
on a given lattice would be the same and equal to that of the corresponding Ising model. Despite the fact that the 
relation is approximate [2 7j . according to available analytical and numerical calculations the critical energy densities 
are indeed, if not equal, very close to each other, whenever a phase transition is known to take place, at least for 
ferromagnetic models on d-dimensional hypercubic lattices. More precisely, the critical energy densities are the same 
and equal to the Ising one for all the 0(n) models with long-range interactions, as shown by the exact solution 0], 
and the same happens for all the 0(n) models on a one-dimensional lattice with nearest-neighbor interactions. As far 
as 0(n) models with nearest-neighbor interactions on cubic lattices with d > 1 are concerned, only numerical results 
are available. According to available data, the transition energy densities are consistent between Ising, 0(2) and 0(3) 
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models on d — 3 cubic lattices, while the critical energies of the ferromagnetic transition of the Ising model and of the 
Berezinskij-Kosterlitz-Thouless (BKT) transition of the XY model on a two-dimensional lattice appear to be only 
slightly different, the difference being about 2% (see Ref. pj and references therein). At present, it is difficult to say 
whether transition energies are really equal for all cases but the BKT one or they are equal only for the long-range 
and the l-d cases and different in all the other cases, with so small a difference that it is not masked by statistical 
errors only in the BKT case. However, the above mentioned results show that assuming the validity of the relation 
as far as the location of the critical energy is concerned — i.e., assuming Eq. ([T]) is a reasonable approximation close to 
the transition energy — gives in general a good prediction for the critical energy itself, and an exact prediction for the 
long-range and the l-d models. It must be stressed that the latter cases are very special, in that the critical energy 
density of the transition equals one of the boundaries of the energy density domain: the lower bound e m i n in the l-d 
case and the upper bound e max in the long-range case. 

Clearly, the fact that a particular prediction made using Eq. (JJJ turns out to be exact does not imply that the 
equation itself is exact. This notwithstanding, it is reasonable to try to understand if in some of the cases where it 
gives the correct prediction for the critical energy the relation (fT]) can be derived with a lesser degree of approximation, 
or even exactly. The aim of the present paper is to show that at least in two particular cases it can be done. The two 
cases that will be considered are the mean- field XY model and the l-d nearest- neighbor XY model, i.e., two n = 2 
representatives of the two classes of 0(n) models where the critical energies are known to be exactly equal to that 
of the corresponding Ising model. We will show that for these two models an expression very similar to Eq. (JlJ, and 
which reduces to Eq. (|TJ) when e — > e c , can be derived exactly in the thermodynamic limit. The technical aspects of 
the derivation strongly rely on the peculiarities of the two models so that we do not see an immediate possibility of 
generalization of the results derived in this paper to generic 0(n) models. This notwithstanding, we are convinced 
that our derivation and results may help in understanding more deeply the relation between 0(n) and Ising models, 
as we shall argue at the end of the paper. 

The paper is organized as follows. In Sec. |TT] the stationary points approach and the approximations introduced 
in [l[ leading to Eq. (fT|) are recalled and discussed. Sees. [Hi] and [TV] are devoted to the explicit derivation of the 
relation between the Ising model density of states and the density of states of the mean-field XY and l-d XY models, 
respectively. In Sec. [V] the results are discussed in a more general perspective, with emphasis on generalization to 
general d dimensional lattices. 



II. STATIONARY POINTS AND DENSITY OF STATES 



Let us now recall the derivation of Eq. (JTJ) made in Ref. As already mentioned in Sec. UJ the approach is an 
"energy landscape" one, i.e., it is based on the study of the stationary points[28| of the Hamiltonian. The importance 
of stationary points in the study of the microcanonical thermodynamics of a system with Hamiltonian H, can be 
understood in an intuitive way as follows. The entropy density s is defined as 

s(e) = ^logu J (e) (2) 

where we have set Boltzmann's constant to unity, and we shall retain this setting throughout the paper. For a system 
with N degrees of freedom described by continuous variables the density of states w can be written as 

where T is the phase space and dT its volume measure, E £ is the hypersurface of constant energy E — Ne, and dH 
stands for the N — I-dimensional Hausdorff measure. The rightmost integral stems from a coarea formula [J]. At a 
stationary point, VH = and the integrand diverges, so that its contribution to uj is clearly important (29j. In the 
following we shall assume that the stationary points contribution to w, or more precisely the contribution of a special 
class of stationary points, is indeed the most important one. 

Let us consider a classical isotropic spin model defined on a lattice (or more generally on a graph) with Hamiltonian 

N N n 

n(n) = - E J ^ - s i = - E J « E s * s ? > ( 4 ) 

— l a—1 

where i and j run over the N lattice sites and the classical spin vectors Si ~ (S},...,S™) have unitary norm, 
i- e -i E"=i (^°) 2 = 1 Vi = 1, . . . , iV. The real matrix dictates the interactions; in case they are long-ranged a 
normalization is understood such as to obtain an extensive energy, using e.g. the Kac prescription [5]. The Hamiltonian 
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(U]) is globally invariant under the 0(n) group. In the special cases n = 1, n = 2, and n — 3, one obtains the Ising, 
XY, and Heisenberg models, respectively. The case n = 1 is even more special because O(l) = Z 2 is a discrete 
symmetry group. In this special case the Hamiltonian (|4| becomes the Ising Hamiltonian 

N 

= - £ JijaiVj , (5) 

ij = l 

where <7j = ±1 Vz. In all the other cases n > 2 the O(n) group is continuous; each spin vector Si lives on an n — 1 
unit sphere 

Let us now consider the stationary configurations of W n ' for n > 2, i.e., the solutions S = (Si, . . . , S at) of the iV 
vector equations 

VH (n) = , (6) 

with the constraint J2a=i( S i) 2 - 

It is impossible to find explicitly all the stationary points (except for some special case, see e.g. @, of for very small 
systems, see e.g. Q and references therein). However, as shown in [![, a particular class of solutions can be found by 
assuming that all the spins are parallel or antiparallel: Sj = ■ ■ ■ = = Mi. In this case the constraints (S*™) 2 = 1 

imply Sf = (Tj Vi, and one finds that the stationary points equations © are satisfied by any of the 2 N possible choices 
of the cr's. The Hamiltonian (j4j becomes the Ising Hamiltonian ([5]) when the spins belong to this class of stationary 
configurations. Therefore we have a one-to-one correspondence between a class of stationary configurations of the 
Hamiltonian (j4j of a 0(n) spin model and all the configurations of the Ising model ([5]), i.e., the Ising model defined on 
the same graph with the same interaction matrix Jij ; the corresponding stationary values are just the energy levels 
of this Ising Hamiltonian. We shall refer to the class of stationary configurations Si = (0, . . . , 0, <jj) Vi = 1, . . . , N as 
"Ising stationary configurations" . There will be also other stationary configurations; nonetheless, the 2 N Ising ones 
are a non-negligible fraction of the whole, especially at large N because the number of stationary points of a generic 
function of N variables is expected to be 0(e N ) [Ell]- 

The above results hold for 0(n) and Ising models defined on any graph. From now on we shall restrict to regular 
rf-dimensional hypercubic lattices and to ferromagnetic interactions > 0. In this case, in the thermodynamic limit 
TV — >• oo the energy density levels of the Ising Hamiltonian ([5]), H^(ai, . . . , <jn)/N Vo"j = ±1, become dense and 
cover the whole energy density range of all the 0(n) models. This fact, together with the above mentioned fact that 
their number is exponentially large in N, suggests that Ising stationary configurations are the most important ones, 
so that we may approximate the density of states a/™- 1 (e) of an 0(n) model in terms of these configurations. To this 
end, let us first rewrite the density of states of our 0(n) spin model as a sum of integrals over a partition of the phase 
space: 



where p runs over the 2 Ising stationary configurations and U p is a neighborhood of the p-th Ising configuration such 
that {Up] 2 p=1 is a proper partition of the configuration space T = (S" -1 )^, that coincides with phase space for spin 
models Since Ising configurations are isolated points in the configuration space of a 0(n) model, such a partition 
always exists. 

Until now we did not introduce any approximation. In Ref. two approximations were introduced to derive Eq. 
(TT]) from Eq. (J7J. Let us review and discuss them, in view of their role in the explicit calculations we shall carry out 
in the next Sections. 

(i) It was assumed that the integrals in Eq. ([7]) depend only on e, i.e., the neighborhoods U can be chosen, or deformed, 

such as 

5{H (n) - Ne)dT = [ 5(H {n) - Ne)dT = g {n \e) (8) 

for any p,q such that H (n) {p) = U {n) {q) = Ne. 

(ii) At a given value of e, the largest contribution to L)( n \e) in Eq. ([7]) is likely to come from those U p such that 

n [n) {p) = Ne, because if H {n) {q) ^ Ne then \\7H<- n) (x)\ f Vx e U q H S e , unless a zero in | VU {n) (x) | comes 
from a stationary configuration which does not belong to the Ising class. Since it was assumed that non-Ising 
stationary configurations could be neglected, only neighborhoods centered around stationary configurations at 
energy density s have been retained in the sum (ff|). 
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As shown in Ref. these two assumptions [3(| immediately lead to Eq. (JTJ) - Both assumptions are needed to arrive 
to Eq. ([T]), and are strictly related to each other. However, these two assumptions might well play a very different 
role. As we shall see in the following sections, in the two analytically tractable special cases, assumption (ii) does not 
hold in general: it holds only when e — > e c . As a consequence, one has to include also stationary configurations with 
energy e' ^ e in the sum. Clearly, if assumption (ii) does not hold, also assumption (i) is of little use as such, since 
also neighborhoods centered around stationary points with energy density different from e have to be included in the 
sum. 

One might then replace assumption (i) with 
(i 1 ) The integrals in Eq. ([7]) depend only on e and on the energy density e' of the stationary point, i.e., 

f 5(H {n) -Ne)dT = G in) (e,e r ), (9) 
Ju p 

for any p such that T-L^ n \p) = Ne' . 
The function would then be related to G^(e,e') by 

g^(s) = G^(e,s). (10) 

Using assumption (i r ) alone, without invoking[31] assumption (ii), one obtains from Eq. j7]) the following expression 
for the density of states of a 0(n) model: 

^(e)=J2^ 1) ^')G (n \e,e'), (11) 

e' 

i.e., a convolution between the Ising density of states uj^ 1 ' and the function G™. Then, in the thermodynamic limit 
N — > oo a saddle-point- like mechanism might single out a value e for e', so that the convolution (fTTj) becomes a 
product: 

wW(e) =w (1) (e)G (n) (e,e), (12) 

where £ is a suitable function of e. If e — e, then using Eq. (fT0|) one recovers Eq. ([T]). This is precisely what happens 
when e — > e c in the two special cases we are going to discuss in the following sections. In sec. [Vj we shall argue about 
the possible generality of this scenario. 



III. THE MEAN-FIELD XY MODEL 



We shall now show that the density of states of the mean- field XY model can be written in the form (fT2|) . with 
e — > e when e — > e c . 

The mean-field XY model Q is a system of N globally coupled planar spins (or alternatively of N globally 
interacting particles constrained on a ring), with Hamiltonian 



1 N 

Hmf = -^ cob^-^) , (13) 



where di € [0, 27r), so that the configuration (or phase) space of the system is the torus T N . This model has a mean- 
field phase transition from a ferromagnetic (or clustered, if one thinks of particles) to a paramagnetic (or uniform) 
phase at e c = e ma x = and has been thoroughly studied being one of the simplest models of systems with long-range 
interactions it belongs to the class (|4]) with n — 2 and Jy ~ 1/N. By introducing the magnetization[32j density 
vector m = (m x ,m y ), where 



1 N 

1 N 



i=l 

(16) 
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we can write the total energy of the system as a function of the modulus m — |m| of the magnetization density: 

Nm 2 



H 



MF 



(17) 



For XY models, Ising stationary points are configurations where the angles #j differ from each other by either or ir. 
Due to the 0(2) invariance of the Hamiltonian, these stationary solutions are not isolated but belong to a manifold. 
We make them isolated by fixing[33| $at = 0, so that the Ising stationary configurations are all the configurations 

$ = where the angles are either or ir, and can be parametrized by the number N„ of angles equal to ir. The 

configurations with given N„ are 



#i = Vi = N n + l,...,N 



(18) 
(19) 



and all the others obtained by permutations of the indices i. The number v(N w ) of such configurations is given by 
the binomial coefficient 



An 



NJ(N-N^\ ' 

while their magnetization and energy density depend only on and are given by 

m(JV w ) = m,(JV T ) = ~* * = 1 - 2n OT , 

m 2 (N„) (N-2N„) 2 (l-2rv) 2 



2JV 2 



(20) 

(21) 
(22) 



where we have introduced the fraction of angles equal to 7r, n.,r = N„/N. 

Given a stationary configuration p — . . . , ^jv}, let us define the neighborhood 



if T?< = 7T 



< 


7T 37T 


2'T 


0< G 


37T 7T 


Y' 2 



(23) 



if i?i = 



so that {Up} 2 =1 is a partition of the phase space T . The density of states lumf of the mean-field XY model can 



thus be written as 



N 



WMF 



where 



G M F(e,iV 7 r) 



3tt/2 



ir/2 



tt/2 



di?i • • • d$ Nw / d&ir w+1 ■■■d$ N 5 [K M f(i?i, . . . ,0*) - We] 



3tt/2 



(24) 



(25) 



We note that v{N v ) given by Eq. (j2"0")) is nothing but the density of states w^f °f the mean-field Ising model 



i N 

H MF ~ ~2jy 2^ 0-1 CT J 



(26) 



as a function of the number of "up" spins a = 1; using the relation (|2"2"|) to obtain the energy density e' of the Ising 
stationary configuration as a function of N w , Eq. (|24[) can be written as 

(27) 



wmf(e) = ^2 w mf( £ ') g mf(e, e') • 



where the sum runs over the energy density levels of the Ising mean- field Hamiltonian (|26l) . so that it is exactly Eq. 
(TTTI) written in the special case of the mean-field XY model. It is important to stress that this result is a consequence 
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of the fact that the energy of a Ising stationary configuration depends only on N v and that all the neighborhoods U p 
with the same N n contribute equally to the sum (|2"4"]l . 

Let us now compute the function Gmf defined in Eq. (|25|) . To make the calculation simpler it is useful to express 
Gmf as a function of m instead of e; one then gets back to e using Eq. (|17p . Since we fixed the magnetization to be 
along the x axis, the function Gmf(w, N„) is given by 



Gmf(iti, N n ) 



3tt/2 



d^.-.d^N, 



tt/2 



r /2 / N \ / N 

/ d& N7r+ x...dD N 6 V] cos tii-Nm] S Vsin^ 



(28) 



Using the integral representation of the Dirac delta distribution, Eq. (|2"5)l becomes 



2 -37T/2 



dtix-d<d Nlt 



exp 



/2 



tt/2 



d'dN^+i-'-d'ff 



N 



iqi I >J cos I?, — Nm 



3tt/2 

exp 



^92 



«?2 ^X] sin ^ ; 



by writing 



^(91,92) 
B(qi,q2) 



p3ir/2 

/ di9 exp [iqi cos + iq2 sin z9] , 

Jtt/2 

pir/2 i-Stt/2 

I d-d exp [iqi cos 1? + z<72 sin$] = / <M exp [iqi cos(d — n) + iqi sm{d — 7r)] 

J3tt/2 Jtt/2 



we get 



G M F(m,iV 7r ) = 



^ \ 2 /»oo /»oo 

27rJ 7 ^ y ^ exp [ N (~ imqi + n 7rlogA(gi,q 2 ) + (1 - n^) log B(q 1: q 2 ))} 



(29) 

(30) 
(31) 

(32) 



The integrals in Eq. ([32]) can be computed with the saddle-point method [10] in the limit N — > 00. The saddle 
point is given by q 2 = e q\ = —27, where 7 £ R satifies the self-consistency equation 



7i( 7 )-L-i(7) 



7i( 7 )+L_i( 7 ) 



w / ( 7 )-Lo( 7 ) +(1 nw) J ( 7 ) + L ( 7 ) 



(33) 



in Eq. (|3"3")1 , 7k (7) are modified Bessel functions of order fc and 7^(7) are modified Struve functions of order k [ill ]. 
We can thus write, in the thermodynamic limit N — > 00, 

G M F(m,n„) = (^j expjiV -m-f + n„ logi( 7 , 0) + (1 - n„) log B( 7 ,0) } , (34) 

where we have written n, instead of N„ since we are in the N — > 00 limit, 7 must be numerically determined solving 
Eq. (|3"3")l , and the functions A and B are 



1(7,0) = tt[/o(7) - £o(7)] , 
5(7,0) = 7r[Zo(7)+io(7)]- 
We can thus write, in the large N limit, the density of states (|24[) as a function of m as 



(35) 
(36) 



wmfH= / efa-Tr exp AT(-m 7 + n,r logj4( 7 , 0) + (1 - n,r) logZ?( 7 , 0) - n,r logn.^ - (1 - n n ) log(l - n n )) , (37) 



where we have neglected the subleading contributions in N. Again, the integral (137)) can be computed with the 
saddle-point method as N — > 00, so that, given m and thus e, only a particular value of (and thus of m' and, in 
turn, of e') is singled out and the density of states wmf assumes the product form (|12p . The particular value of n w 
which is singled out is the one such that the exponent in Eq. (|37p is maximum; it has to be computed numerically. 
The saddle point on Eq. (|37[) singles out a value m of the magnetization such that 



u>MF(m) =u/ 1 )(m) Gmf(to,to) 



(38) 
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h(m) 




FIG. 1: Numerical results for the function h(m) defined in Eq. (|39jl for the mean-field XY model. The red (dotted) part of the 
curve is obtained by interpolation (see text). 

In order to show that the value of m as a function of m converges to m as m — > m c , where m c = is the critical value 
of the magnetization, in Fig. Q] we plot the function 

h[m) = m — rh . (39) 

Figure Q] shows that h — > as m —> 0, so that the density of states umf(m) is such that 

wmf(w) -> w (1) (m) g M F(m) , (40) 

where <?mf(w) = Gmf(wi> m), for m — > m c . More precisely, h appears to be a linear function of m as m -) 0, 
h(m) oc —m. When m — > 1 the numerical procedure we used to compute h(m) had some convergence problems. 
Since m — 1 implies h(m) = and n T = 1, to avoid these numerical problems the curve plotted in Fig. Q]in the 
range m £ [0.97, 1] has been evaluated interpolating the numerical results obtained for m < 0.97 with the constraint 
h(l) — 0. The interpolating curve is drawn in red and in dotted style in Figure[TJ We stress that the part of the curve 
relevant to the phase transition is that in the opposite limit, m — > 0, where the numerical procedure easily converges. 
We can now go back to the energy, using s = — m 2 /2, and write 

u MF {e)=u J W{e)GMF(e,i), (41) 

where e — > e as e — » e c — 0. One can thus write, as e — > e c , 

wmf(e) -> (e) 0mf(e) , (42) 

where 3mf(^) = Gmf(s,£), for e — > e c . Figure shows the function 

V (s)^h(V^2~£) = £-e; (43) 

as e — > e c = 0, 77(e) cx —\J~^e. Since \rj(e)\ is the difference between the energy e singled out by the saddle point 
and the energy e at which the density of states is calculated, it somehow measures also the "distance" between the 
function Gmf(e, e) and the function gyiF (^O — Gmf(e, £)■ From Fig.[5]we see that this difference reaches its maximum 
(roughly equal to 1.2 x 10~ 2 ) around the center of the energy density range. Comparing this value to the width of 
the energy range itself we see that this difference is at most of the order of 2%. 

IV. THE ONE-DIMENSIONAL XY MODEL 

Let us now consider the one-dimensional XY model, which is a system of N planar spins with nearest-neighbor 
coupling, described by the Hamiltonian 

N-l 

Ui d = -J2 cos - **) . ( 44 ) 

i=l 
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FIG. 2: Numerical results for the function 77(e) = h(— \/2e) denned in Eq. (|43|) for the mean-field XY model. The red (dotted) 
part is obtained by interpolation (see Fig. [T] and text). 



where, as in the mean-field XY model, -di <E [0, 2tt), so that the configuration (or phase) space of the system is the 
torus T N . This model does not have a bulk broken symmetry phase; it is ordered only in its state of minimum energy. 
Hence, the phase transition from a ferromagnetic to a paramagnetic phase occurs at e c = e m i n = —1, and at zero 
temperature. It belongs to the class ^ with n — 2 and Jij = 1 for i and j nearest-neighbors and zero otherwise. 
As we shall see in the following, also for this model the density of states can be written as 

w w (e) = «W(e)Gi <I (e,e) ) (45) 

where, in this case, a/ 1 - 1 is the density of states of the one-dimensional Ising model 

iV-l 

n il = - E > ( 46 ) 

i=l 

and e — > £ as e — > e c — e m - m . One can thus write, as e — > e c , 

w ld (e)->wW(e)sid(e), (47) 

where gid(e) — G\d(e,e), for e — > e c . The derivation follows very closely that of the mean-field model, with a few 
differences that will be underlined. 

Let us fix i? jv =0, and leave open the boundary condition at the other side of the chain. As in the mean- field case, 
the Ising stationary configurations are those where the angles § are either or n. However, their energy is no longer 
parametrized by N v . On an Ising stationary configuration, the energy can be written as 

Hu@i, ■ ■ .Jn-i) = = 2N d - N + 1 , (48) 

where is the number of the domain walls in the configuration, i.e., the number of flips between 7^ = and $ — ir 
(and viceversa) along the chain. This implies that one can no longer use the definition (|^|) of the neighborhoods 
Up to build the partition of the configuration space, because this would imply that stationary points with the same 
energy would give different contributions. 

Let us then change variables from ($1, . . . , t?at) to {x\, . . . , x^) as follows: 

x k = tf fe+ i-tf fc iffc = l,...,JV-l, 

(49) 

XN = $N = . 

In the new variables the Ising stationary points are still such that Xk. = or x^ = n, but now the energy is given in 
terms of the number of re's equal to 7T, because the number of domain walls Nd is precisely that number. One can 
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thus define the partition of the configuration space using the neighborhoods U p defined as 

if Xi = 7T 



Xi e 



7T 37T 

2'T 



37T 7T 

~2~ 5 2 



(50) 



if X; = 



and write the density of states of the 1-d XY model as 

N-l 



uu(e) = 51 v{N d )G ld (e,N d ) 
N d =a 



where 



(N-l)\ 



N d \(N-N d -l)\ 



(51) 



(52) 



is the number of Ising configurations with N d domain walls, i.e., the density of states a/ 1 -*^') of the one-dimensional 
Ising model with energy density 



2N d - N + 1 
N 



and 



Gid(e, N d ) 



r-Zix/l rir/2 I N-l 

/ dxx ■ ■ ■ dx Nd I dx Nd+ i ■ ■ ■ dx N -i 6 - V] 

Jtt/2 J3tt/2 \ k=1 



cos a;^ — iVe 



(53) 



(54) 



The computation then proceeds following very closely what already done for the mean-field case. The 1-d case is even 
simpler, because one can directly compute G\ d as a function of the energy density, without the need to consider it as 
a function of the magnetization. Using the integral representation of the 6 and integrating on the x variables we can 
write in the large N limit 



1 f°° 

G ld (e,n d ) = — / dqexp{N[-iqe + n d logb(q) + (1 - n d )\oga(q)}} , 

^ J-oo 



where n d — N d /N and the functions a and b are given by 

r/2 



a(q) = / dxexp (— iq cos x) , 

J3n/2 



b(q) 



r/ 
r 3Tr/2 



I dx exp (—iq cos x) 

Jtt/2 



Performing again a saddle point with q = —ij we get, in the N —¥ oo limit, 

1 



G ld (e,n d ) = -expjiV -je + n d log 6(7) + (1 - n d ) log 0(7) | 



where 



a( 7 ) = A( 7 ,0), 
6(7) = 5(7,0), 

with .A and B given by Eqs. (|35|) and (|36|) . respectively, and where 7 satisfies the self-consistency equation 

n , Ji(7) - ^-1(7) , /i(7)+£-i(7) 

E = (1 - rid) r ; s r , s + n d ■ 



4(7) - ^0(7) 



/o(7)+io(7) 



(55) 

(56) 
(57) 

(58) 



(59) 
(60) 



(61) 
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We can thus realize that Eqs. (l58l) and (|6Tj) coincide with the same equations derived for the mean-field case, i.e., 
Eqs. (J31J and fl33J), provided 

( m ^ 1 C 62 ) 

The latter reflect the fact that in the l-d case the transition occurs at the minimum value of e instead of at the 
maximum. 

From now on, the calculation of u>id(s) is exactly the same as that of wmf(w), with the substitutions (|62[) . A given 
value of fid of rid will be singled out, which corresponds to an energy density e via Eq. (|53[) . We thus obtain 

w u (e)=«W(e)Gi <I (e,e) ) (63) 

where e — > e as e — > e c = e m - m ; more precisely, defining the function 

C(e) = e - e = fc(m = £ + 1) , (64) 

where /i(m) is the function (f3T)| defined for the mean-field model, we have that £ — > when e — > e c = e m i n = — 1, 
and in particular £ oc — (1 + e) for e close to e c = — 1. If one plots C as a function of e one thus obtains exactly 
the same curve reported in Fig. [TJ with the horizontal axis shifted so that e € [—1,0]. Since |/i(m)| is maximum for 
m ~ 0.75, the function \((s)\ reaches its maximum value (roughly equal to 0.15) around e ~ —0.25; the maximum 
difference between e and e in this case is around 15% of the full energy density range, larger than in the mean-field 
case. 



V. CONCLUDING REMARKS 

The present paper has been mainly devoted to discuss the validity of the relation |T]). put forward in [l|, in the 
special cases of the mean-field and l-d XY models. We have shown that a slightly more general formula, Eq. (|12p . 
holds, which reduces to the previous one in the limit e —> e c . 

The two models we have dealt with are very special and both of them are exactly solvable in the microcanonical 
ensemble. This feature is crucial for the derivation we have presented. As a consequence, a generalization of the 
computations to 0(n) models with short-ranged interactions on a d-dimensional lattice is not straightforward at all, 
the difficulties being similar to exactly solving their thermodynamics in the microcanonical ensemble. 

This notwithstanding, we can learn something from these results. The present work confirms that, as already noted 
in [l[, Eq. |T]) can not be exact for a generic value of the energy density; at most, it could be valid for e = e c . Indeed, 
it was already pointed out in [l[ that the relation ([1) could not be exact for a generic 0(n) model, since the specific 
heat critical exponent a of a 0(n) model would then have the correct sign, but the wrong absolute value. More 
precisely, Eq. {]]) implies that if ai is the microcanonical specific heat exponent of the Ising model on a given lattice, 
then the microcanonical specific heat exponent of the 0(n) model on the same lattice and with the same interactions 
is a — —a j, regardless of n. In d = 3, for instance, this yields the correct sign of the 0(n) exponents, because ai > 
so that a < 0; the 0(n) specific heat is not divergent, but cuspy at the transition. However, the absolute value of 
the exponent is wrong, because it should depend on n, as shown by well-established results for the 0(n) universality 
classes [e|. It is worth noting that, here and in the following, we are dealing with the specific heat critical exponents 
defined in the microcanonical ensemble: these are related to the usual critical exponents a defined in the canonical 
ensemble by a = a/(l — a) [HI, so that microcanonical results can be easily carried over to the canonical ensemble 34]. 

The result a = — aj follows from Eq. |T]) by assuming that the function g( n \e) is a generic function which does 
not contain any explicit information on the phase transition, i.e., is analytic with a generic Taylor expansion. If we 
proceed in an analogous way assuming that Eq. (TT2)) holds for a generic 0(n) model, we still find the correct sign of 
the specific heat critical exponents as with Eq. [I] but we do no longer have any contradiction with the known results 
on the values of the exponents. Indeed, assuming that G(") (x, y) is a generic (i.e., analytic) function because it should 
not contain any information about the phase transition, it can be shown that the critical exponent a of the continuous 
model can be any real number in [—1,0). This range of values is in agreement with known results [l2[; moreover, 
although it does not predict a precise value of a, it still correctly implies that the specific heat of 0(n) lattice spin 
models does not diverge for n > 1. The details about the predictions of Eqs. |T]) and (|12[) as to the critical exponent 
a are reported in Appendix [A] The latter observations, together with the observed near-equality of critical energies 
of 0(n) models defined on the same lattice at different n [l|, suggest that the relation (fT2j) might have a more general 
validity than being restricted to the two special cases considered in this paper. 

Although the above argument shows that a more general validity of Eq. (|12p as an exact result can no longer 
be excluded on the basis of the predictions for the exponent a, it would imply, as in the case of Eq. ([1} that was 
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already discussed in [l|, a "Patrascioiu-Seiler" scenario in d — 2, i.e., the presence of some kind of phase transition 
in two-dimensional 0(n) models also for n > 2 [l4|, EH- This scenario, althou gh n ot ruled out by rigorous results, is 
believed to be unlikely on the basis of numerical simulations (see e.g. Ref. [I(ljl8l|). 
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Appendix A: Specific heat critical exponent 

In Sec. El we have discussed the implications of Eqs. (TTJ) and (fT2"j) in case they would exactly hold. Here we give the 
details about the predictions on the specific heat critical exponent a obtained by assuming that the density of states 
has the form given by Eqs. ([1} or (fT2")l . respectively. Let us recall that, in the microcanonical ensemble, the specific 
heat is defined as 

<*.>~Jg$, (AD 

where s(e) is the entropy density and the temperature is defined as T(e) = l/s'(e). With s'(e) and s"(e) we denote 
the first and second derivative of the function s(e). 

Let us consider a short-range 0(n) model and assume the relation (JTJ) holds as an equality. We assume in the 
following that the phase transition occurs for a value of the energy density in the interior of the domain of the entropy 
density [33. Without loss of generality, let us shift the energy density e such that e c — 0. The entropy density of the 
continuous model can then be written as: 

s(e) =*/(e) + log /(e), (A2) 

where here and in the following we use the notation s/(e) instead of s^'(e) for the entropy density of the Ising model, 
to avoid possible misunderstanding with derivatives. We also omit the symbol (n) indicating which 0(n) model we 
are considering because our arguments do not depend on it. Finally, we denoted gWfe) 1 /* by /(e). 

Let us now consider, for the moment, only energy densities larger than the critical one, i.e., e > 0. Three facts are 
relevant for the following: 

1. we consider < ai < 1, i.e., the case d > 2. Moreover, because the critical temperature of the Ising models is 
finite, Sj(e) cx e a ' for e-)0 + . 

2. s'(e) is finite around e = because the critical temperature of the continuous model does not vanish at the 
transition. 

3. we assume /(e) is analytical, consistently with the discussion in Sec. |TJ We can then expand /(e) in a Taylor 
series around e = 0. 



Inserting Eq. (|A2[) into Eq. (jAll) . we get 



2 



C ( £ ) = ~ —^—72 • (A3) 



„ll( \ , 3"(e) 
S A £ ) + /( E ) 



m 



Using the expansions described above around e = 0, neglecting the higher order terms and expanding the fraction, 
we obtain 

C(e) ~ a+ + b+ e°" (e -> 0+), (A4) 

where a+ and 6+ are constants whose exact value is irrelevant to our purposes. We can repeat the same calculations 
for e < 0, obtaining the same result as in Eq. (|A4[) but for that e — > —e and that the constants may be different. 
Hence the specific heat close to e = is 

C{e)~a±+b±\e\ ai . (A5) 
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We then obtain the result stated in Sec.|V] the specific heat of the continuous model does not diverge at the transition 
and the critical exponent a of the continuous model is related to the one of the Ising model via a — —aj. 

With a similar reasoning we can also deal with the case in which we consider Eq. (fT2j) to be exact. As before, we 
start by considering e > 0. Assuming Eq. (fl"2"j) holds as an equality, the entropy density of the continuous model is 

s(e) = Sl (e)+f(e,e(E)), (A6) 

where we denoted by f(e, e(e)) the function (1/N) log g( n ' (e, e(e)) . In this case, / is a function of two variables: again, 
we assume it is analytic and expand it around e — 0, such that 

f(x, y) ~ / + fix + fry + f 3 xy + f 4 x 2 + f 5 y 2 + f 6 x 2 y + f 7 xy 2 + f s x 3 + f 9 y 3 , (A7) 

where x and y are shorthands for e and e and the f^s are constants whose exact value is irrelevant to our purposes. 
At variance with the previous case, e(e) contains some information about the transition because it vanishes for e — > 0; 
we should then admit the possibility of a singular dependence on e, writing e(e) oc e e with 9 > for e — > + . 
Using the information on the behavior of s'/(e) around e — and integrating two times, we get 

s/(x) ~ a + aix + a 2 x ai+2 , (A8) 

where the a^'s are suitable constants. Inserting Eqs. (IA7|) and (IA8|) into the equation for the entropy of the continuous 
model, Eq. (|A6|) . we get: 

s(e) ca Q + ai e e + a 2 e d ^ + f + he + f 2 e 9 + f 3 e e+1 + f,e 2 + f 5 e 29 + f 6 e 2+e + f 7 e 1+29 + / 8 e 3 + f 9 e 3e . (A9) 
Taking the first and the second derivative of the previous expression and renaming the constants, we obtain 

s'(e) ~ he 9 - 1 + b 2 e e[ai+2) - 1 + hi + h 2 e e - x + h 3 e e + h 4 e + h^ 2 ^ 1 + h 6 e e+1 + h 7 e 26 + h 8 e 2 + hge 36 ' 1 , (A10) 
and 

s"(e) ~ cis 6 ~ 2 + m 3 £ 0_1 + m 4 + m 5 e 28 ~ 2 + m 6 s e + rnje 20 ^ 1 + m 8 e + m 9 e 38 ~ 2 . (All) 

The quantity is unknown. However, since the specific heat of the continuous model does not vanish at the transition, 
the above expressions imply the constraint 9 > 2. Moreover, if 9 > 3, the linear term in Eq. (|A11|) would dominate. 
Hence the range of values for 9 to be considered is 9 £ (2, 3]; if 9 > 3 or 6 = 2, the leading behavior of s"(e) would 
be the same as that given by Eq. (|All[) with 9 = 3. 

The leading behavior of Eqs. (IA10[) and (IA11[) is then s'(e) ~ hi + h^e and s"(e) ~ 7714 + c\E 6 ~ 2 . Inserting these 
results into the expression (IA1[) for the specific heat, we obtain 

C{£) ^_{hi+!H^_ + + d+£ e-2 (£ ^ 0+) _ (A12) 

Repeating the same calculations for e < and combining the result with Eq. (|A12[) we obtain the behavior of the 
specific heat close to the transition, 

(7(e) ~ c± +d± \e\ e ~ 2 . (A13) 

The above expression, together with the above bounds on 9, shows that the specific heat of the continuous model 
does not diverge and its critical exponent a is determined by 9, which is model dependent. Varying 9 in its allowed 
range we obtain a € [—1,0). 
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